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PREFACE 


The work described in this report was performed by the Space Sciences 
Division of the Jet Propulsion Laboratory. 
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ABS t RACT 


Skylab-4 X-5 thermal data of The Ueysers area waa analysed to determine 
the feasibility of using midday Skylat images to detect geothermal areas. 

The hottest ground areas Indicated on the Skylab image corresponded to south- 
facing barren or sparsely vegetated slopes. Thermal well D 4, a geothermal area 
approximately 15 by 30 m coincided with one of the hottest areac indicated by 
Skylab. However, this area could not be unambiguously distinguished from 
the other areas which are believed to be hotter than their surroundings as a 
result of their topography, and mlcrometeorologlcal conditions. A simple 
modification of Watson's thermal model was performed and the predicted tempera- 
tures for the hottest slopes using representative values was in general agreement 
with the oh' -ved data. It is concluded that data from a single midday Skylab 
pass cannot be used to locate geothermal areas. 
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INTRODUCTION 


Ceotheraal exploration has received considerable attention within the past 
several years as a result of the Increased der nd for energy. Areas favorable 
for such exploration have, in general, been recognized from field reconnaissance 
and by remote sensing from aircraft. Although remote sensing techniques 
have been utilized for almost a decade, until the advent of Skylab, no non- 
military satellite had sufficiently high thermal or spatial resolution to detect 
these areas. 

Theoretically, the ideal times for detection of such areas are postdawn 
and at sunset (Watson. 1975). The effects of thermal inertia of surface mat- 
erials are minimized at these times. Unfortunately, Skylab data acquisition 
times do not coincide with these optimal times. The present study was under- 
taken to determine if data from a Skylab near-noon pass could be used to detect 
a known geothermal area. 

The study area chosen was the Ceysers geothermal field in Sonoma County, 
California, which is one of the few producing geothermal fields in the world. 
This area has been studied extensively using both ground and a*** . survey 

techniques, (see Allen and Day, 1927; Cal. Dlv. Mines and CeoJ . . ’66; McNltt, 

1968; Moxham, 1969; Stanley, et al., 1973; U.S.C.S., 1973; Chapman, 1975), 
and more recently by the Martin-Marietta Company (1974) using Skylab nighttlre 
data. 

Data analyzed was acquired by one experiment carried out in the Skylab 
missions, the S 192 multn'oectral scanner. This scanner Imaged the ground 
in 13 wavelength bands in the visible, reflective infrared and thermal Infrared. 
Data was recorded in a 68.5 km swath about the spacecraft ground track. 

Two different detectors were used for the thermal infrared .hannel (10.2 
to 12.5 um) at various times during the 3 missions. The Y-3 detector was used 
for the first two missions, and part of the third. It had a theoretical thermal 
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resolution of 0.7*C, but due to excessive noise, this resolution was not realised, 
and the detector was replaced by the X-5 detector during the third ■lsslon. The 
actual NEAT of this detector was approxlnately 0.67*C. Only data from the X-5 
detector was used in this study. 

DESCRIPTION OF THE SITE 

Geography and Topography 

The Ceysers geothermal atea la located in Sonoma County, California, 100 km 
north-northwest of the city of San Francisco. An area approximately 15 by 12 km 
was chosen for study (Figure 1). 

The topography of this region is characterized by northwest trending valleys 
and ridges. Relief varies from 400 to 1100 meters (Figure 1). Drainage for the 
area is through Big Sulfur Cre?k, Squaw Creek, and Little Sulfur Creek, which all 
feed Into the Russian River to the southwest. The climate Is Mediterranean, and 
most c f the precipitation occurs in the winter; summers are hot and dry. Vege- 
tation consists of seasonal grasses, scrub bushes and some trees. 

The main geothermal activity occurs in the valley occupied by Big Sulfur 
Creek and along Geyser Canyon, trending north from The Ceysers Resort. The south- 
facing slopes along Big Sulfur Creek are barren or sparsely grass covered with 
occasional scrub brush and trees along drainages, whereas north-facing slopes 
are predominantly tree covered. 

Geology 

The Geysers area is underlain by the o irass lc-Cretaceous Franciscan Assemblage, 
a complex eugeosyncllnal sequence of chert, greenstone, graywacke and sedimentary 
breccia, and sedimentary and tectonic melange. The Franciscan Assemblage is 
overlain by Jurassic (?) metamorphosed ultramaflc rocks and serpentine. 

Extrusive volcanic rocks of Fliocene-Plelstocene age occur near the study 
area. These include the Clear Lake series, which extends to within 8 km of The 
Geysers area, and rhyolite flows which cap Cobb Mountain just northeast of the area. 
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The area has undergone complex faulting, trending predominantly northwest, 
producing a serlea of horats and graben. The Individual horcts and graben are 
2-3 km wide and extend for up to 15 km (Figure 2). 

The Geysers la one of numeroua thermal areas within one of these northwest- 
trending graben. McNltt'a (1968) Interpretation Indicates A km of vertical dis- 
placement along this graben. It Is possible that this displacement has provided 
channels for the ascension of heated fluids, which could be the source for geo- 
thermal activity in the area. 

Recent detailed geologic mapping of The Geysers region by McLaughlin (1974) 
shows that the area Is extensively covered by Quaternary landslide deposits. 
Hydrothermal alteration is prevalent, predominantly along fault zones, and 
transcends lithologic boundaries. 

Geothermal Development of Region 

Geothermal activity was first discovered in The Geysers area in 1847, and 
the region soon became a well-known health resort. The f 4 rst wells were drilled 
in 1921 for the purpose of generating electric power. By 1925, eight wells were 
completed; however, because of the lack of a market for steam-generated electricity, 
the project was discontinued. 

Development resumed in 1958 when the Pacific Gas and Flectric Company con- 
tracted with Magma Power Company and its partner Thermal Power Company to build 
a steam-electric power plant. This plant (units 1 and 2 In Figure 3) went Into 
final production In 1963, and now generates 24,000 kilowatts of electricity by 
steam. 

To date, a total of six power plant units have been built, five are opera- 
tional. Each unit consists of a generating plant approximately 20 by 40 met»rs, 
and a cooling unit, approximately 40 by 100 meters. Steam Is supplied to these 
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planta from numerous nearby • •ells (Figure 4). Total power p aductlon In 1974 
was approximately 400,000 kilowatt of electricity, enough to supply the needs 
of a city half the sixe of San Francisco. 


DIGITAL DATA PROCESSING 

Skylab digital data of The Geysers geothermal area was obtained from 
Johnson Space Center on computer compatible magnetic tapes. Reflected and em- 
itted radiation was recorded by the S-192 multlspectral scanner in thirteen 
discrete wavelength band.*, of which six were used (Table 1). The Instantaneous 
field of view of one picture element (pixel) is approximately 75 by 75 meters. 
These date were processed by the Image Processing Laboratory at JPL. 

TABLE 1. Wavelength of Selected Skylab S-192 Channels 


SD0 1 

Wavelength (pm) 

03 

.56 - .61 

07 

.68 - .78 

09 

.78 - .88 

11 

1.55 - 1.75 

15 

10.20 - 12.50 

19 

.98 - 1.03 

21 

10.20 - 12.50 


A geometric correction was applied to the digital data, unmodified in the de- 
lineated study area, to compensate for skew caused by the earth's rotation during 
acquisition of the image. The image was then resampled digitally :o Improve the 
registration to established topographic base maps of The Geysers region. Contrast 
in the rectified image was then Increased to distribute measured brightness levels 
(DN) over the entire dynamic range of the photographic film used for display of 
the digital image. This process increases the investigator's ability to discern 
contrast variations in the picture. Prior to contrast enhancement a probability 


4 


JPL Technical Memorandum 33-728 


density function (PDF) vaa created which deacrlbea the frequency of occurrence 
of DN within an image. An appropriate contraat atretch waa derived baaed upon 
the PDF. The PDF waa alao uaed to id** tlfy the brlghteat 100 and 300 pixela, 
corresponding, in the thenaal image, to the hotteat polnta (Flgurea 5 and 6). 

Halloing of two reflective channela waa uaed to reduce the effect of topog- 
raphy on raolometry, thus enhancing apectral Information. Ratlolng conalata 
of dividing, pixel by pixel, the DN valuea of the lmagea, and encoding theae 
ratioa aa DN to form a new image. 

A detailed deacrlptlon of theae technlquea la given by BUllngaley and 
Goett (1973). 

DATA ANALYSIS 

Preliminary Analyala 

Aa a first step in the recognition of the hotteat ground areaa within The 
Geysers region, a probability density function (aa described above) of ground 
temperatures was produced from the Skylab-4 X-5 near-noon pass. No anomalous 
hot spots were detected from an analysis of the thermal channel data (SDO-21). 
Rather, ground temperatures were normally distributed. If geothermal areas 
are to be distinguishable, their temperature must be significantly higher than 
their surroundings. Therefore, computer enhancement techniques, as discussed 
previously, were used to delineate the hottest 100 and 300 ground areas (pixels). 
These results are displayed In Figures 3 and 6. The hottest areas are black, 
decreasing ground temperatures are In diminishing shades of gray, and the coolest 
areas are white. Although there appears to be an extreme contrast between the 
hot areas and their surroundings, this is due to the enhancement procedure; actual 
temperatures differ less than 1°C between the black and the darkest gray regions. 
Stream channels and other pronounced physiographic features can be recognlred 
and matched to the same features in the topographic map (Figure 1). 
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Proa Inspection of Figures 5 end 6, It Is appsrent thst sll of the hottest 
sress occur on south-fsclng slopes. It Is also apparent that the cluaters of 
the top 300 pixels are primarily an expansion of the clusters of the top 100, 
Indicating that these are probably true hot areas and not noise In the data. 
Inspection of aerial photographs Indicates that almost all of these areas occur 
on bare ground. These observations would be expected as a reault of normal 
midday aolsr heating. 

Comparison of the top 100 and 300 pixels with the most recent detailed 
geologic map of The Geysers region (Figure 2) indicates tnat the hottest areaH 
appear to be Independent of 1 thology. Comparison of these areas with published 
data on the location of known geothermal areas and related man-made structures 
(Figures 3 and 4) shows only slight coincidence. 

The obvious relationship between topography and solar heating Introduces 
the otrong possibility that all the hot areas were those most favorably oriented 
to recelv'' solar radiation. Therefore scatter diagrams of the temperature 
of the top 300 pixels (expressed in DN) vs. angle of incidence of solar radia- 
tion on these slopes were produced. The scatter of points wa*» too large to 
determine a relationship between these variables. Separate scatter diagrams 
were produced for the major geologic units in order to determine if this scatter 
was due to the variability in surface properties of different lithologies. No 
relationship could be established. Despite these negative results, a relation- 
ship should exist between ground temperature and angle of Incidence of solar 
radiation. The most probable cause for this scatter was that an Insufficient 
range of temperature data was being analyzed. Therefore, a more thorough 
quantitative analysis, with o larger data base, was undertaken. 

Quantitative Analysis 

As a first step In the quantitative analysis, an attempt wa» made to super- 
impose precisely the Skylao image on topographic maps (The Geysers and Asti, 


6 


JPL Technical Memorandum 33-728 


California 7 1/2' quadrangles) containing the atudy area. Due to the lack of 
reaolutlon of proalnent topographic featurea within the Skylab data and minor 
distortions, Individual plxela could not be located exactly. Becsuse of thla 
uncertainty, a 3 x 3 pixel average was uaed to reduce thla aource of error. The 
3x3 pixel grid waa placed on the topographic maps and the alope angle and 
azlauth were measured In two aelected reglona for 924 9-pixel groupa. Theae 
two reglona contained the majority of the hlgheat temperatures and the known 
geothermal areaa. 

To minimize the effects of vegetation, only the 145 9-pixel groupa were 
analyzed which coincided with barren or sparsely grass covered areas aa Identi- 
fied on aerial photographs. The cosine of the angle between the normal to these 
slopes and incident solar radiation was plotted against the average DN value 
(Figure 7). The expected relationship between temperature and solar radiation 
is apparent, In contrast to the previous analysis of only the top 300 pixels. If 
areas with significant geothermal heating are present, their distribution of 
points on Figure 7 should lie above and to the left of the main distribution. 

The large scatter of data precludes unambiguous Identification of any geothermal ly 
heated area in the study region. Furthermore, refinement by analyzing each geo- 
logic unit separately did not reduce the scatter. 

A contingency table was constructed and a chi square test was performed to 
determine the association between the bare areas and the clusters of the pixels 
representing the 100 hottest temperatures which were contained in the two 
selected regions (see Table 2). The observed association between bare areas 
and hot spots is indicated in Table 2a. If these variables are Independent, 
the expected relationship would be that given In Table 2b. A chi square test 
was performed on the tabulated values to determine the significance of the ob- 
served relationship. Chi square is calculated from 

X* ‘ l (0-E)» . 

i - 1 E 


JPL Technical Memorandum 33-728 


7 


where 


0 ■ Observed 
E ■ Expected. 

For the observed sssoclstlon X 2 “ 102.6. According to statistical tables, 
(Arkln and Colton, 1963), the expected X* value for 1 degree of freedom at 
Pm is 10.8. This implies that the observed association between bare areas 
and hot spots would occur by chance much less than once in a thousand. 


(a) Observed 


T able 2. Contin g ency Table 


Cround Cover 



Bare 

Not Bare 

Total 

Hot 

24. 

7. 

31. 

Not Hot 

169. 

1143. 

1312. 

Total 

193. 

1150. 

1343. 

(b) Expected If 

Independent 

Cround Cover 



Bare 

Not Bare 

Total 

Hot 

4.45 

26.55 

31. 

Not Hot 

188.55 

1123.45 

1312. 

Total 

193. 

1150. 

1343. 


The association between bare areas and hot spots, and the importance of 
the topographic control on solar heating required the introduction of computer 
data handling techniques to evaluate and then remove the effects of these vari- 
ables. Digital topographic data of The Oeysera region were provided by Honeywell 
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Marini? Systems Division, West Covina, Calitornla. These data were extracted 
from the 1:250,000 topographic map of the Santa Roaa quadrangle, and consisted 
of the ground elevation in feet on a 208.33 foot (63.5 meters) orthogonal grid. 
For this established grid, the slope angle and azimuth were computed from the 
digital data. The amount of solar radiation incident o.. these slopes was com- 
puted from an existing atmospheric model (Cates et al., 1971). The results 
are displayed in Figure 8. The intensity of incident radiation is displayed 
in varlc is gray levels with white being the most Intense. The Skylab data 
were converted to the topographic coordinate system to compare the thermal 
data with the solar irradiation data. It is apparent that the slope informa- 
tion derived from the digital topographic data lacks the resolution necessrry 
to match accurately the detail provided from Skylab. 

Because the effects of vegetation on surface heating have not been quanti- 
fied techniques refined by Rowan et al. (1974) were used in an attempt to 
identify vegetated areas and remove them from the data set. Ratio pictures 
of the varloun Skylab spectral bands were produced for this purpose. Figure 9 
is a ratio of SDO 3 (.56 to .61 um) to SDO 7 (.68 to .76 pm) and is typical of 
the ratio pictures produced from the Skylab data. Unfortunately, this procedure 
also er'iances the noise and it is apparent that in this data the njlse obscures 
any new meaningful information. 

Despite the lack of success in removing completely the effects of topog- 
r-i'uy ...id vegetation, the analysis provided valuable insight into the magnitude 
ji these effects. 

DETECTABILITY OF GEOTHERMAL AREAS BY SATELLITE 

Ground temperatures are determined by measuring the flux ground infrared 
raoiation which reaches the spacecraft sensors. The sensor Integrates the flux 
over a resolution element and provides a signal corresponding to an area weighted 
average temperature after calibration. For the purposes of this investigation. 
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a simple calculation was performed to determine the site and temperature 
of a surface geothermal feature which could be detectable by Skylab. The 
effects of atmospheric absorption were Ignored. 

The flux emitted by a surface of temperature T is EOT 1 ' , where C la the 
emlsslvlty of the surface, which Is assumed here to be 1.0, and 0 Is the 
Stefan-Boltzmann constant. 

The Skylab thermal imaging system yields an average temperature over 
a ground area, A, of approximately 75 meters by 75 meters. If a smaller 
area, a, at a hotter temperature, T^, exists within the resolution element 
A, the measured value (T g ) will be the total flux from the various temper- 
ature sources 


A 0 T 4 - a O T 4 + (A - a) 0 T 4 

s H g 


where T^ Is the average temperature of the remainder of the pixel. Solving 
for a v have 


a 


A(T 4 
s 

(T 4 
v H 




This equation can be used to establish the size and temperature of an 

anomalously hot area which could be detected by Skylab sensors. This 

anomalously hot area must be hot enough and/or large enough to raise 

the integrated temperature sensed by the satellite, T g , significantly higher 

than the temperature of the surrounding pixels. A background temperature, 

T , of 280°K is assumed for the normal temperature of the ser.sed pixel and 

g 

its surroundings. The size and temperature relationship for detectability 

at various levels of T above the background temperature, T , are shown in 

s g 

Figure 10. For example, a hot area 16 meters by 16 meters at a temperature 

of 360°K would raise the sen&wJ '■emperature , T , by 5° to 285°K. Clearly, 

s 
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geothermal areas must he hot and/or large to be unambiguously detected by 
Skylab sensors. Published reports and maps of The Geysers area Indicate no 
naturally occurring hot springs or steam vents that fulfill the detectability 
requirements. There is, however, one large anomalously hot region, which is 
detectable: thermal well #4 (see discussion on Field Reconnaissance). 

In contrast to the convective tranrport of hot fluids, which creates only 
isolated small hot spring and steam vents in The Geysers area, surface areas 
considerably larger than the areal extent of a pixel can be heated by con- 
duction. Although the satellite will be able to detect these elevated tem- 
peratures as long as they are above the noise level of the sensors (approximately 
1°C), these areas cannot be positively identified as geothermal regions without 
consideration of all the other components of the surface heat balance equation, 
These include solar heating, terrestrial radiation, atmospheric radiation, 
latent heat of evaporation, and sensible heat transport into the atmosphere. 

All of these are one to several orders of raagn .tude larger than the geo- 
thermal heat flux which will be induced by the known temperature gradients 
in The Geysers area. With a temperature gradient of 20*C/100 m and a con- 
ductivity of .01 cal/°K cm sec the geothermal heat flux would be 20 heat 
flow units (HFU). Hase (1971) in studies at Yellowstone National Park, Wyoming, 
has suggested that a geoti'emal heat flow must be greater than 240 HFU to 
be detectable by any method. Based on a theoretical model Watson (1975) has 
shown that with only a single temperature measurement the error in the heat 
flux can easily vary from 250 to 1000 HFU. Therefore, even though the possi- 
bility exists of enhancement of surface temperatures in The Geysers area due to 
geothermal heating by conduction, recognition of such areas is highly unlikely. 


1 1 


JPL Technical Memorandum 33-728 


FIELD RECONNAISSANCE 


Field reconnaissance was undertaken In The Ceysers regie i to try to deter- 
mine the nature of the heating processes occurring at the location of the 
pixels representing the hottest temperatures. One cluster of these pixels 
(cluster A In Figure 5) coincided with the highest thermal anomaly seen by 
Moxham (1969) in daytime and nighttime infrared Images of The Ceysers region. 

This area includes a large geothermal region of boiling surface waters and 
steaming ground about a blown-out steam well, thermal well #4. The highest 
temperature value within the cluster occurs at the well location. 

The hottest ground temperatures within this region occur in a hydrothermally 
altered area approximately 15 by 30 meters. Although measurement" of subsur- 
face temperatures at a depth of a few centimeters v«re consistently near boiling, 
surface occurrences of boiling water and steam are patchy. Between these very 
hot areas, suri^ce temperatures were only slightly to moderately elevated due 
to apparent evaporative cooling. 

Several drilled vents in the immedi te area are forcefully venting steam, 
which combined with the steaming ground, produced a large steam cloud over 
the entire area, to a height of about 100 meters. Temperature of the venting 
steam at the smallest drilled hole was 125°C at the top of the well casing, 
decreasing to 70°C two meters above the well casing. 

Although steam wells, small steam vents and minor surface springs were 
observed in the immediate surrounding areas corresponding to the remaining 
pixels in cluster A, it was impossible to determine if there was sufficient 
geothermal heat flow to enhance measurably the surface temperature. These 
regions are steep, barren or sparsely grass covered, have south-facing slopes, 
and receive substantial solar heating. 
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Based upon these field observations, it appears certain that the geo- 
thermal heati. _ about thermal well #4 was detected by Skylab. The actual 
temperature measured by Skylab was 25°C for the pixel coinciding with thermal 
well #4 t which is 4°C hotter than its surroundings. This temperature is rea- 
sonable in view of the effects of evaporative cooling, masking of the ground 
by steam clouds, and averaging of temperatures over an area larger than the 
anomalously hot region about thetaal well #4. 

An extensive network of steam pines leads from the wells to power plants 
and cooling units in many parts of the study area. None of these structures 
are hot enough or large enough to be identified by Skylab. 

The hottest temperature within The Geysers region did not occur near 
thermal well #4 nor near any other known geothermal site. It is located 
above Squaw Creek (Cluster B in Figure 5) on a steep, barren, southeast- 
facing slope comprised of dark greenstone. Between outcrops, slopes were 
thinly covered with dry grass. No surface indications of geothermal activity 
were found. 

Clusters C and D are located above Squaw Creek and also contained some of 
the hottest temperatures (Figure 5). They were similar in all respects to 
Cluster B. 

Cluster E is located on the northwest side of CobL Creek (Figure 5). The 
southeast-facing slopes are more heavily grass covered, less steep, and have 
fewer rock outcrops than the Squaw Creek locations. The slopes are concave 
and form a sheltered valley. Therefore, local meteorological effects may be 
significant in the heating of this region. Minor occurrences of warm springs 
indicate geothermal activity in the general vicinity, which may also enhance 
the heating. 
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Cluster F on Figure 5 Is located at the abandoned Cloverdale Mine. Although 
it does not Include one of the highest 10 temperatures, this cluster was inves- 
tigated because of the mapped occurrence of hydrothermally altered rocks, 

(Figure 2) and reports of gas. The slopes are ste^p, southeast facing 
and sparsely grass covered. Extensive mine dumps of chert mantle the slopes 
beluw the Cloverdale Mine. The hottest temperature in the cluster coincided 
with graywacke outcrops above the dumps. No geothermal activity was noted. 

All *he aforementioned locations except thermal well 14 had concave, 
sheltered slopes, though to a lesser degree than at Cobb Creek. 

Areas of known but unpublished geothermal activity were investigated and 
no correlation could be found with the top 300 pixels on the Skylab image. 

All of these areas were either heavily vegetated or located on flat, west- 
facing or north-facing slopes. 


THERMAL MODELING 

Thermal modeling was undertaken to explain the occurrence and 
distribution of the hottest ground areas within The Geysers region. This 
investigation stemmed from the observation that the hottest ground tempera- 
tures were, without exception, on south-facing slopes. In addition, the 
effects of topography on solar heating in an area of rugged terrain, such as 
The Geysers, should be expected to overwhelm the very subtle effects of geo- 
thermal heating of the surface. 

The thermal model used in this study was developed by Watson (1971; 
Watson, et al., 1971; Watsot., 1973, 1975). This model explains the surface 
temperature behavior of the gre ind through a diurnal cycle. Watson, (1975) 
indicated that a sampling frequency of three or four measurements during a 
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single diurnal cycle would be required to unambiguously Identify geothermal 
anomalies below 100 heat flow units under Ideal conditions. However, If 
data can be acquired only once during a diurnal cycle, the optimum time 
would d*. either at sunset or shortly after sunrise. Neither of these 
time constraints was met with the Skylab-4 overflight data of The Geysers 
region. Nevertheless, an attempt was made to explain the distribution 
of hot rpots within The Geysers region by using representative values for the 
necessary parameters in Watson's model. 

The temperature values and locations of the pixels w^re transferred 
to a topographic map of The Geysers and the local dips and azimuths were 
read for the 15 hottest slopes in the ccene.* 

An examination of the S-190 B photograph of The Geysers area showed 
that each of the slopes had a high albedo and was barren of any signifi- 
cant amount of vegetation. The values for the 15 slopes (dip, azimuth, 
and temperature averaged over the slope) were Incorporated into the Watson 
model in an attempt to find a single set of values of thermal inertia 
and albedo which would most closely match the observed temperature of 

these slopes in the S-192 data. The parameters for the initial model 

-2 -U 

were as follows: thermal inertia • .035 cal cm sec , albedo • .55, 
emissivlty ■ .95, sky temperature ■ 260°K, cloud cover factor ■ .20. These 
parameters are consistent with a mixture of soil and rock with a high 

albedo. Thermal inertias of soils generally vary from .010 to .040 cal 

-2 - 1 * -2 
cm sec , rock from .040 - .090 cal cm sec . 

*The selection of slopes was based on the presence of a cluster of several 
pixels whose individual temperatures were at least 18°C. The geographic 
limits of the pixels accepted for any given slope were determined by any 
significant change in dip or azimuth of the slope (in each case, between 
9 and 15 pixels). 
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Although the temperature extremes tor the lb slopes were only b.¥ s C 
for the Skylab data, the model predicted a temperature difference of 10.0*C. 

Any variation In the parameters which would reduce the predicted tempera- 
ture difference necessitated the use of unrealistic values of cloud cover, 
sky temperature, r emloslvlty. A plot of azimuth versus midday tempera- 
ture generated In the Initial model for various dips Is shown In Figure 11. 

A modification of the Watson model was made which assumed that slopes 
across the valleys from the measured slopes occulted a significant fraction 
of the sky. Because north-facing slopes were nearly at sky temperature, 
only a small modification of sky temperature was necessary for south-facing 
slopes. However, the south-facing slopes were found to be 30 to 35°C hotter 
than the sky and a large adjustment In sky temperature was necessary for 
north-facing slopes. East-and west-facing slopes were intermediate In 
value. 

A simplistic modification of the initial model was assumed: that each 
of the 15 slopes was looking at a 30 degree opposite-facing slope which 
occulted one-third of the sky, and sky temperatures were adjusted accord- 
ingly. No additional correction was made for reradlation and reflection between 
surfaces, although the effects were clearly important. 

Figure 12 is a plot of temperature generated In the modified model 
versus azimuth for various dips. The lowest temperatures attained in this 
model are -2.5*C, a value which Is compatible with the S-192 temperatures. 

The data from this model was compared to temperatures from the Skylab 
thermal channel. Figure 13 is a plot of the angle between the sun and 
surface normal, versus temperature for the 15 hottest slopes. A 1.1 # C 
correction factor was added for atmospheric attenuation. There is very 
close agreement between observed and predicted temperatures. Consequently, 
explanation of the hot areas does not require geothermal heating. 
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ANALYSIS OF NORTH-FACINC SLOPES 


The highest temperatures in The Ceysers region coincide with bare, 
south-facing slopes. To reduce the effect of solar radiation, analysis 
of nighttime Images is desirable. Data for the only nighttime pass of 
The Ceysers region was examined, but because of the exceptionally high 
noise level, was not considered. 

Therefore, an attempt was made to make use of shadowed steep-slope 
areas as simulated nighttime Images. The temperatures of these north- 
facing slopes were extracted for the daytime X-5 thermal channel (SDO 21). 
The highest 100 temperatures within these regions were set to black and 
displayed in picture form (Figure 14). Within the accuracy of pixel 
location, the ground areas with highest temperatures coincided with mar- 
ginally sunlit ridges or slopes. No correlation was found between these 
pixels and geothermal areas on the north-facing slopes. 
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SUMMARY AND CONCLUSIONS 


Skylab X-5 thermal data for a daytime pasa over The Geyaera geothermal 
field waa analyzed to determine the feaalblllty of using Skylab to detect 
geothermal areaa. The hotteat ground areaa were Identified on the Skylab 
Image and were located on topographic maps, aerial photographs, and In 
the field. 

In conclusion, the following statements can be made: 

1) The hotteat areaa found on the S-192 thermal lmagea in The 
Ceyaera region are all found on aouth-faclng a 1 ope a which are 
barren or aparaely vegetated. 

2) No statiatlcally algnlficant association exlata between the 
hotteat areas on Skylab thermal Images and known geothermal 
areas in The Ceyaera region. This result la consistent with 
the temperature and spatial resolution capability of the S-192 
scanner, and published field data of the area. 

3) Although thermal well #4 coincides with a pixel which represents 
hot ground, it cannot be unambiguously distinguished from areas 
believed to be hotter than their surroundings as a result of 
their topography, micrometeor logical conditions, and insolation. 

4) A simple modification can be made in the Watson model which 
considers reradlatlon from opposite facing slopes. This mod- 
flcatlon produces a close agreement between observed and 
predicted temperatures without requiring geothermal heating. 
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ng. 1 . Topographic map of a portion of the Kelseyviile 15' Quadrangle showing The 
Geysers geothermal region 







area after MrLaughlin (1974) 


2. Generalized geologic map 


‘ * : 



Fig. 3. Location map of geothermal wells and power plants in The Geysers area (after State of 
California Division of Oil and Gas, 1974) 
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Fig. 5. The Geysers region at 10. 2—12. 5 urn, computer enhanced to show the pixels 
representing the 100 hottest temperatures 
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Fig. 6. The Geysers region at 10. Z— IZ. 5 pm f 
representing the 300 hottest temperatures 


computer enchanced to show the pixels 


rv 


» 



cos z* 


Fig. 7. Graph of DN versus cosine of the angle between 
slope normal and sun (Z ) 
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Fig. 9. Ratio picture of SOO 3 (. 56 to . 61 |im| to SDO 7 (. 68 to 
. 78 >irr>) 



Fig. 10. Detectability of thermal anomalies as a function of 
size and temperature 
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PARAMETERS 


Azimuth (*) 

Fig. 11. Temperature versus azimuth for various dips (initial model). 
Note that the lowest temperature attainable is -13*C, a value which is 
considerably colder than the lowest temperature measured from the 
S- 1 92 data 
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Azimuth (•) 


Fig. 12. Temperature versus azimuth for various dips (modified 
model). Note that the lowest temperature attainable is -2. 5 * C, a 
value which is compatible with the S-192 data 
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Angle between sun ond surface normal 

C) 


Fig. 13. Temperature versus angle between the sun and 
surface normal for the 15 slopes used in the model 
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Fig. 14. Extracted north-lacing slopes of The Geysers region. 

are black 


The 100 hottest ground areas 



